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Abstract 

The asymmetries of both light-flavor antiquark d{x) — u{x) and strange-antistrange s(x) — s(x) 
distributions of the nucleon sea are considered with more details in the effective chiral quark 
model. We find that the asymmetric distribution of light-flavor antiquarks d(x) — u{x) matches the 
experiment data well and that the asymmetry of strange and antistrange distributions can bring 
about 60-100% correction to the NuTeV anomaly of sin which are three standard deviations 
from the world average value measured in other electroweak processes. The results on the correction 
to the NuTeV anomaly are insensitive to the inputs of the constituent quark distributions and the 
cut-off parameters. The ratios of d(x)/u(x) and s(x)/~s(x) are also discussed, and it is found that 
the ratio s(x)/~s(x) is compatible with the available experiments with an additional symmetric sea 
contribution being considered effectively. 
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I. INTRODUCTION 



The physics community has been puzzled since the NuTeV Collaboration reported the 
anomaly that the value: sin 2 6 W = 0.2277±0.0013 (stat) ±0.0009 (syst), which was measured 
in deep inelastic scattering (DIS) of (anti-)neutrino from the nuclear target with possible un- 
certainties and errors having been considered [lj], is deviated from the standard model (SM). 
The Weinberg angle (or weak angle) 9 W is one of the important quantities in the electroweak 
theory, and can be determined from various experimental methods, such as atomic parity vi- 
olation, W and Z masses, elastic and inelastic neutrino scattering, and so on. At present, the 
world accepted value is sin 2 6 W = 0.2227 ± 0.0004. The NuTeV Collaboration measured the 
value of sin 2 6 W by using the ratio of (anti-)neutrino neutral-current to charged- current cross 
sections on iron jl|, which is closely related to the Paschos-Wolfenstein (P-W) relation Q|: 

°cc °cc z 

The validity of this relation is based on the fundamental assumptions of isoscalar target, 
charge symmetry and strange symmetry s(x) = s(x), which give rise to a large number of 
theoretical and experimental investigations on this issue. Possible sources of the NuTeV 
anomaly beyond the SM have been discussed jjj. However, before speculating on the pos- 
sible new physics, one should first check carefully the standard effect and the theoretical 
uncertainties within quantum chromodynamics (QCD). 

Here, we make a broad review of the main mechanisms which have been used to investigate 
the NuTeV anomaly. The P-W relation holds for a nuclear target provided that the nucleus is 
in the isoscalar state, which means that various strong interaction effects must cancel out in 
the ratio. But the actual targets used in the neutrino experiments are usually non-isoscalar 
nuclei with a significant neutron excess, such as the iron target in the NuTeV experiment. 



Kumano 



investigated a conventional explanation in terms of a nuclear modification caused 
by the difference between the u- and d-valence distributions when the charge and baryon- 
number conservation for nuclei are considered. He and his collaborators estimated this effect 
to the NuTeV sin 2 6 W value by using a \ 2 analysis method to reproduce nuclear data on the 
structure function F2 and the cross-section of Drell-Yan processes, and noticed that the 
effect is not large enough to explain the whole NuTeV anomaly in their later work And 
Kulagin |(| not only took into account the neutron excess correction to the P-W relation 
for a neutron-rich target, but also discussed other corrections to the P-W relation caused 



by other nuclear effects: the Fermi motion, nuclear binding, nuclear shadowing, etc. In 
addition, nuclear effects were completely estimated and the shift to the NuTeV anomaly 
was given with some uncertainties in Ref. []| by using a particular nuclear x-rescaling model 
to describe the structure functions in (anti-)neutrino-nucleus DIS. Also, there are other 

n 

suggestions |8j from a conservative point of views. 

Another fundamental assumption is charge symmetry which is a more restricted form 
of isospin invariance involving a rotation of 180° about the "2" axis in isospin space, or 
more specifically, the isospin symmetry for the u <-» d exchanges between protons and 
neutrons However, there have been some discussions about the correction to the P- 

n 

W relation due to this symmetry breaking. Sather [H| first pointed out that the charge 
symmetry violation (CSV) should largely affect the extraction of sin 2 6 W from the neutrino 
collision and gave an estimation about it. The similar result was given in Ref. that 
CSV effect should reduce about one-third of the discrepancy between the NuTeV result and 
the accepted average value of sin 2 ^. On the contrary, Cao and Signal calculated the 
nonperturbative effect of CSV within the meson cloud model framework and showed no 
contribution to the NuTeV anomaly. At this stage, we still have no direct experimental 
evidences to point to a substantial violation of CS in parton distributions. Only there are 
some up limits in experiment for the CSV, because the effect of CSV usually confuses with 
the asymmetry of s and s distributions. 

The asymmetry of strange and antistrange quark distributions may imply an additional 
source for the anomaly. The topic of strange content in the nucleon sea is one of the most 
challenging issues in the hadron physics, especially for its connection to the proton spin 
problem and to the NuTeV anomaly. For the part of nucleon sea generated through gluon 
splitting g — ■> qq in perturbative QCD, CP symmetry is expected in the quark and antiquark 
distributions in the leading order. It has been argued recently that there is a strange- 
antistrange asymmetry in perturbative QCD at three-loop jlflj . however, the perturbative 
QCD alone definitely predicts a non-vanishing and Q-dependent value of strange-antistrange 
asymmetry which would have a very small contribution to the extraction of sin 2 ^. So 
the cause of a sizable asymmetry of the nucleon strange sea should be of nonperturbative 
origin jl^l I15I Q|. Furthermore, it is highly probable that the momentum distributions 
of strange and antistrange are different in the nucleon sea, although the total number of 
strange quark and antiquark occurring as quantum fluctuations must be precisely equal 
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to conserve the strangeness quantum number. The consequence on the NuTeV anomaly 
from possible asymmetry of strange and antistrange quark distributions in the nucleon sea 
was examined by Cao and Signal utilizing the meson cloud model 15, Q] and the 



result is fairly small and has no significant effect on the NuTeV result. Oppositely, it was 
shown that the light-cone baryon-meson fluctuation model proposed by Brodsky and 
Ma Q| can describe the second moment S~ = x[s(x)—s(x)}dx and produce an asymmetry 
of s(x) and s(x) distributions which could remove roughly 30-80% of the discrepancy between 
the NuTeV result and other determination of sin 2 ^. Also, in Ref. [3], the asymmetric 
s — s distribution treated in a different framework of non-perturbative hadronic K + A 
fluctuations is found to reduce the NuTeV result to only about 2 standard deviations from 
SM, while the realistic behavior of the strange sea distribution [xs(x) + xs(x)]/2 is also 
reasonably reproduced. Kretzer et al. obtained QCD correction by calculating the 
next-to-leading-order (NLO) neutrino cross sections and studied the shift of sin 2 9 W which 
is closely correlated with s(x) ^ s(x) parton distribution function sets and isospin violation 
with some uncertainties. It is also noticed that the contribution to the NuTeV anomaly 
changes sign from "— " to "+" when S~ is calculated [sjj by using the Lagrange Multiplier 
method. A more reliable theoretical analysis in the chiral quark model has been give in 



Ref. 



22| , in which it is shown that the NuTeV anomaly can be accounted for by at least 60% 



without sensitivity to the inputs of constituent quark distributions and cut-off parameters. 
Recently, Wakamatsu j3] gave also a theoretical analysis based on the flavor SU(3) chiral 
quark soliton model by introducing a parameter of the effective mass difference between 
strange and nonstrange quarks and predicted a fairly large asymmetry of s- and s-quark 
distributions which would solely resolve the NuTeV anomaly, in similar to the conclusion in 



Ref. [22j. In addition, Szczurek et al. |24| discussed more earlier the dressing of constituent 
quark with a pseudoscalar meson cloud within the effective chiral quark model by including 
the effect of SU(3) f symmetry violation explicitly. They found more pions and kaons in the 
nucleon than what predicted by traditional meson cloud model and pointed out that the 
effective chiral quark model may lead to a sizeable asymmetric strangeness content of the 
nucleon sea. 

As we know that the nucleon sea has received attention for a long time because of the 
abundant phenomena which are away from the naive theoretical expectations. For example, 
the Gottfried sum rule (GSR) 2^] violation is related to the light-flavor antiquark asymmetry 
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in the nucleoli sea and has been measured in several experiments. At the same time, there 
are some models, such as: the meson cloud model |26j, |27] in which the relative success 



was obtained in both the difference and the ratio of d and u by including the perturbative 



component; Pauli blocking effect which first suggested by Field and Feynman 



28| : and 



several different types of chiral models and so on, proposed to explain the violation of 
Gottfried sum rule. And another case is the proton spin crisis j^J and its connection to the 



strange content in the nucleon sea [30] . Usually, the distributions of s and s are assumed to 
be symmetric, but this is neith er p roved theoretically nor experimentally. Moreover, some 



teitn er p rov 

0,0fl 



measurements and analyses |31l.l32ll33l] show that the distributions of strange and antistrange 
quark may be asymmetric. A joint fit [3| to the neutrino charged-current cross sections 
and charged lepton structure function data was made by the CDHS Collaboration and some 
improvements were achieved in the fits if asymmetry in the strange sea was allowed. 

In this paper, we present the more detailed calculation of the asymmetric distributions 
for both the light-flavor antiquarks d(x) and u(x) and the strange content in the nucleon sea 
by using the effective chiral quark model along with our previous work [2^, and find that 
the distributions of d(x) — u(x), d(x) + u(x), x(s(x) +s(x)) and s(x)/s(x) match well with 
the experimental data, when additional symmetric sea contributions being considered effec- 
tively by taking into account the difference between model results and data parametrization. 
Furthermore, it is should be pointed out that the asymmetry of s(x) and s(x) distributions 
could remove roughly the NuTeV anomaly by at least 60%, and it is more remarkable that 
this conclusion is insensitive to the different inputs of the effective chiral quark model. 



II. THE SEA CONTENT IN THE EFFECTIVE CHIRAL QUARK MODEL 



The effective chba! qu a r k mo de,, established by Weinberg ffl, and deve!oped by Mancta 
and Georgi [37| . has been widely accepted by the hadron physics society as an effective theory 
of QCD at low energy scale. The effective chiral quark model has an apt description of its 
important degrees of freedom in terms of quarks, gluons and Goldstone (GS) bosons at 
momentum scales relating to hadron structure. There has been a prevailing impression that 
the effective chiral quark model is successful in explaining the violation of GSR (which was 
first detected by the New Muon Collaboration at CERN by the analysis of Eichten, 
Hinchliffe and Quigg B| ■ Also, this model plays an important role in explaining the proton 
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22j shows that the strange- 



spin crisis |29| by Cheng and Li (^J. A recent study by us 
antistrange asymmetry within the effective chiral quark model could explain the NuTeV 
anomaly also. It is the purpose of this paper to provide a more detailed analysis on the 
strange-antistrange asymmetry and its connection to the light-flavor sea asymmetry in the 
effective chiral quark model. 

The chiral symmetry at high energy scale and its breaking at low energy scale are the 
basic properties of QCD. Because the effect of the internal gluons is small in the effective 
chiral quark model at low energy scale, the gluonic degrees of freedom are negligible when 
comparing to those of the GS bosons and quarks. In this picture, the valence quarks which 
are contained in the nucleon fluctuate into quarks plus GS bosons, which spontaneously 
break chiral symmetry, and any low energy hadron properties should include this symmetry 
violation. The effective interaction Lagrangian is 



(2) 



where 



d 



(3) 



is the quark field and = + igG^ is the gauge- covariant derivative of QCD, with 
standing for the gluon field, g standing for the strong coupling constant and qa for the 
axial-vector coupling constant determined from the axial charge of the nucleon. and 
are the vector and axial-vector currents which are defined by the following forms: 




(4) 



where £ = exp(zll//), and II has the form: 
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(5) 



With the expansion for and in powers of IT//, it gives = + 0(II//) 2 and A^ = 
idfjl/ f + 0(11/ f) 2 , where the pseudoscalar decay constant is / ~ 93 MeV. So the effective 
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interaction Lagrangian between GS bosons and quarks in the leading order becomes 

Ln q = - 9 -j^d^y^. (6) 

The framework that we use in this paper is based on the time-ordered perturbative theory 
in the infinite momentum frame (IMF). Because all particles are on- mass-shell in this frame 
and the factorization of the subprocess is automatic, we neglect all possible off-mass-shell 
corrections. In this framework, we can express the quark distributions inside a nucleon as a 
convolution of a constituent quark distribution in a nucleon and the structure functions of a 
constituent quark. The light-front Fock decompositions of constituent quark wave functions 
have the following forms 

\U) = \/Z\u ) + a n \dTT + ) + ^\u7T°) + a K \sK + ) + -^=|w/>, (7) 

v2 V6 

\D) = y/Z\d ) + a>0 + ^\dn°) + a K \sK°) + ^\d V ), (8) 

V 2 V6 

and the corresponding picture are shown in Fig. ^ Here, Z is the renormalization constant 
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FIG. 1: The dressing of the constituent U and D quarks with Goldstone bosons. 



for the bare constituent quarks which are massive and denoted by \uq) and \do), and |a a | 2 
are the probabilities to find GS bosons in the dressed constituent quark states \U) for an up 
quark and \D) for a down quark, where a — 7T, K, n. In the effective chiral quark model, the 
fluctuation of a bare constituent quark into a GS boson and a recoil bare constituent quark 



can be given 



4l|: 



<?j0) = / —Pja/i(y)q%{-)- 
Jo y y 



(9) 



This process also can be described in Fig. Efa). In Eq. Q, Pj a /i(y) is the splitting func- 
tion which gives the probability for finding a constituent quark j carrying the light-cone 
momentum fraction y together with a spectator GS boson a and having the following form: 



9Am^ 2 f 2 (rrij - rriiy) 2 + 



57T 



f ' J V(i-y)K-^J a 



where m^m^m a are the masses of the i, j-constituent quarks and the pseudoscalar meson 
a, respectively, m = (m^ + rrij)/2 is the average mass of the constituent quarks, and M? Q is 
the invariant mass square of the final states: 



Mf a = ^l + 4±^. (11) 

y i-y 

In this paper, we also adopt the definition of the moment of splitting function: 




FIG. 2: (a) A constituent quark qi fluctuates into a Goldstone boson a plus a recoil constituent 
quark qj. (b) A Goldstone boson emits a qi and a qj . 
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(12) 



with the first moment (Pj a /i) — (Paj/i) = (Pa) = |«o| 2 Q| It is conventional that an 
exponential cutoff is used in IMF calculations because the integral in Eq. requires a 
momentum cutoff function at the quark-GS boson vertex. Commonly, 



9 A -> #A ex P 



4A 2 



(13) 



with g' A = 1 following the large N c argument |42|, however, g' A = 0.75 was taken in the origi- 
nal work 3?J • And A is the cutoff parameter, which is determined by the experimental data 
of the Gottfried sum and the constituent quark mass inputs for the pion. The experimental 



S 



value of Gottfried sum is 0.235 ±0.0026 |38(, which is different from the value 1/3 predicted 
by the naive parton model in which both the light-flavor sea and the isospin between proton 
and neutron are SU(2) symmetric 3|. Employing the quark distributions of the effective 
chiral quark model, we get the Gottfried sum determined by the difference between proton 
and neutron structure functions: 

- 1 dx, 



S, 



Gottfried 



[Fi{x) - F2{x)\ 

o x 
1 f 1 dx 



L u(x) + u(x) — d(x) — d(x)] 
3 Jo x 

\{Z - \ (Pit) + (P K ) + l <P„» 

1(1- 2 (P n )). 



(14) 



From above equation, we gain the values of Gottfried sum given in table H] corresponding to 
different A,,-, from which we can find that the appropriate value for A^ is 1500 MeV, when 
the value of Gottfried sum matches well with the experimental data. At the same time, 
(P w ) = 0.149, (P K ) = 0.085, (P v ) = 0.063 and Z = 0.682. But for K and 77 mesons, the terms 
(Pk) an d (P v ) in the Gottfried sum cancel out those terms in Z = 1 — j(P n ) — (Pk) — !(-F??)> 
so the value of A can not be determined from Eq. (fUjl or experimental data. It is natural to 
assume that the cutoffs are same for ir, K and 77 mesons in the effective chiral quark model: 



A^ = A, 



1500 MeV 



4l| . which is different from the traditional meson cloud model. 



TABLE I: The values for A and Gottfried sum 
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When probing the interna 
process in the following form 



structure of GS bosons in Fig. Efb), we can also write the 



4i| 



dy!dy 2 ,x y x 

Vk/a{ — )Paj/i{ — )qi{y2), 

yi V2 yi 2/2 



(15) 



where P a j/i(x) = Pj a /i(\ — x), Vk/ a (x) is the quark k distribution function in a and satisfies 
the normalization Vk/ a (x)dx = 1. Because the mass of rj is so high and the coefficient 
is so small so that the fluctuation of it is suppressed, the contribution is not considered in 



our calculation in this paper. From Eqs. J7J) and JHJ), we can have the quark distribution 
functions of the nucleon by using the splitting function Eq. (jlOj) and the constituent quark 
distributions u and d which are normalized to 1, 

U(x) = Zu (x) + P U7T -/ d 0d O + K/7T+ ® Pn+d/u ® Mo + ^PunO/u ® «0 
+ ® Pr:+ S /u ® U + ^K/ttC ® (Pr<V« ® M + PnOd/d ® do), 

d(ar) = Zd (a:) + P d7r +/ U ® « + Vd /vr - ® P w - u /d ®d + -P d7r o /d ® d 

+ Vrf/jro <g> P*o a /d (g d + jVd/wO ® (PrO u/u ® u + P^ (8) d )- (16) 
Here, we define the notation for the convolution integral: 

P®q= f^-P{y)q{-l (17) 

Jx y y 

and 

WP«,= f* /^^Pd-^fe). (18) 

Jx Vi J yi V2 yi y2 

In the same way, we can have the light-flavor antiquark and strange quark and antiquark 
distributions: 

u(x) = V~ /7r - ® P^-u/d ® d + -Va/ w o ® (P^ou/u ® u + P,rO d/(i ® d ), 

d(^) = V 1/7r+ <g> Pr+d/^ ® MO + J V d/TT° ® (PrO«/« ® M + PrOrf/d ® <fo) , 

s(ar) = P S ic+/u<8>Mo + P s ico/<i<8)do, 

s(x) = Vs/K+ ® PftT+s/u ® n + (g) P^s/d g) d , (19) 

where 



2Vu/w° — 2Vu/tt° = 2V rf / 7r U = 21^0 



~K(x,/4lo), (20) 



and 



Vj/K+ = Vs/K°, Vu/K+ = Vd/K°- 

From above equations, we can reexamine the valence quark distributions u v (x) = u(x) —u(x) 
and d v (x) = d(x) — d(x) which satisfy the correction normalization with the renormaliza- 
tion constant Z. The parameterizations of parton distributions for mesons are taken from 
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GRS98 M and the sea content in the nucleon is not considered here, 

K(z,/4lo) = 1.052x- a495 (l + 0.357v^)(l-a;) - 365 , (21) 

V u/K+ (x,ft 2 NLO ) = V d/K o(x,ft 2 NLO ) = 0.540(1 - x) ai7 K(x,/4 LO ), (22) 

Vs/ K +(x, /iNLo) = Vs/k°(x, /4lo) = /4lo) ~ K*/Jir+(z, A*nlo) • ( 23 ) 

Thus, we can calculate the asymmetries of light-flavor antiquark distributions d{x) — u(x) 
and strange-antistrange quark distributions S~ = x[s(x) — s(x)]dx within the effective 
chiral quark model. In this present work, we choose m u = = 330 MeV, m s = 480 MeV, 
m n ± = m n o = 140 MeV and m K + = m K o = 495 MeV. Constituent quark (CQ) model 
distributions [44 1 and CTEQ6 parametrization Q| as two different kinds of constituent 
quark distributions of inputs are adopted. The CQ model distributions have the following 
forms: 

"°W = B[c 1 + 1 ,c 2 1 + c 2 + 2l - C, ' 1 -^ +a+ ' 

B[c 2+ l|2c 1 + 21 ^ (1 -^'- (24) 



where B[i, j] is the Euler beta function and c\ = 0.65 and C2 = 0.35 adopted from Ref. |44 14( 
are independent of Q 2 with the number and momentum sum rules: 



1 f i 

u (x)dx = 2, / d (x)dx 
Jo 

1 r l 



/ xuo(x)dx + / xdo(x)dx = 1. (25) 
Jo Jo 



It is pointed out that there are other different values for c\ and C2 suggested by Ref. 
And the CTEQ6 parameterizations are: 

u (x) = 1.7199x- a4474 (l -x) 2 ' 9009 exp[-2.3502x](l + exp[1.6123]x) 1 - 5917 , 

d (x) = 1.4473x-°- 3840 (l - x) 4 ' 9670 exp[-0.8408x](l + exp[0.4031]x) 3 - 0000 . (26) 

As is well-known, the quark distributions measured by experiments at certain Q 2 include 
not only non-perturbative intrinsic sea but also perturbative extrinsic sea Q| . The intrinsic 
sea quarks are relatively Q 2 independently (or less dependently) multi-connected with the 
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valence quarks in the nucleon, but the extrinsic ones are generated mainly from the gluon 
splitting. Because the correlations between the intrinsic and extrinsic sea quarks are less 
Q 2 dependent, we give the distribution of d(x) — u(x) at low energy compared with the 
experimental ones at high scale shown in Fig. El From the figure, we find that the calcu- 
lated results using the effective chiral quark model match the experiments well with two 
different sets (CQ model and CTEQ parametrization) of inputs for the constituent quark 
distributions. The difference between the results from two different inputs may reflect some 
uncertainties due to the evolution effect which should not cause a big influence on the con- 
clusion. Besides that, we also calculate the strange sea distributions of x(s(x) +s(x)) and 
x(s(x) — s(x)) within this mechanism. For the distributions of x(s(x) + s(x)), we give the 
comparison of the results from the effective chiral quark model calculations (with two dif- 
ferent sets of inputs) and the NuTeV data parametrization j^] obtained from the LO fits 
to the dimuon production cross section in z^Fe and V^Fe deep inelastic scattering in Fig. 0] 
From the figure, we can find that the results of the effective chiral quark model calculations 
are lower than the parametrization of NuTeV data at arbitrary x, which may be caused by 
the non- considered symmetric strange sea content which will be discussed in Sec. IIVI 

In the effective chiral quark model, it should not be the same for the momentum of 
strange quark (s) coming from the recoil constituent quark and that of anti-strange quark 
(s) originating from the virtual GS bosons K. Using Eqs. (|16[) and (|19p. we can obtain the 
distributions for x5 s (x) shown in Fig.EJwith S s (x) = s(x) —s(x). From the figure, we can see 
that the magnitudes of x5 s (x) with CQ model as input are almost twice larger than those 
with CTEQ parametrization as input, however, the contributions to the NuTeV anomaly, 
which will be shown in the next section, are similar. 

III. APPLICATION TO THE NUTEV ANOMALY 

In 1973, Paschos and Wolfenstein suggested to consider the observable Eq. by 
measuring the ratio of neutral-current to charged-current cross sections for neutrino and 
antineutrino on isoscalar target to reduce the uncertainties from the productions of the heavy 
flavor quarks. In Eq. a u ^ c (erj^) is the integral of neutral-current inclusive differential 
cross section for neutrino (antineutrino) over x and y, and the same to Gqq ((Jq 1 ^). The 
NuTeV Collaboration was inspired by the P-W relation to measure the Weinberg angle 
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FIG. 3: Distributions for d(x) — u{x) for A„- = 1500 MeV. The solid and dashed curves are the 
model calculation results in the effective chiral quark model with the constituent quark (CQ) model 
and the CTEQ6 parametrization as inputs for uq(x) and do(x), respectively. The data are from 
HERMES (Q 2 = 2.3 GeV 2 ) Q and E866/NuSea (Q 2 = 54 GeV 2 ) experiments 
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FIG. 4: Distributions for x(s(x) + s(x)). The solid (dotted) curves are the calculated results in 
the effective chiral quark model with the CTEQ6 parametrization (CQ model) as input. The thick 
solid (dotted) curves are for Ak = 900 MeV and the thin ones are for Ax = 1100 MeV, respectively. 
The dashed curve is for the parametrization of NuTeV data 3] • 
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FIG. 5: Distributions of x5 s {x) with 6 s (x)=s(x) — s(x) in the effective chiral quark model, when 
the inputs for uq(x) and do(x) are both constituent quark (CQ) model (thick curves) and CTEQ6 
parametrization (thin curves) with Ak = 900 MeV (solid curves) and 1100 MeV (dotted curves ). 



sin 2 6 W using high statistics separated neutrino and antineutrino beams which arise from 7r 
and K mesons decay following the interaction of 800 GeV protons at FNAL. Because the 
charged- current events contain a final muon that penetrates substantially farther than the 
hadron shower, neutral- and charged- current events are distinguished by the event length 
in the counter. The resulting interaction events must be deposited between 20 GeV and 
180 GeV in the calorimeter and will be found in the NuTeV detector. The most general form 
of the differential cross section for neutral-current interactions initialed by (anti-)neutrino 

n 

on nucleon target is |51j : 



2v(p) 



dz 
a 



NC 



dxdy 



ITS 



a 



2 sin 2 9 W cos 2 6 W M Z 



Ml 



M 2 Z + Q? 



■[xyF?{x,Q 7 



+(1-2/ 



xym 2 N 



V 



Ff(x,Q 2 )±(y-^)xF*(x,Q% 



(27) 



And similarly, for the (anti-)neutrino-nucleon charged-current reaction, the differential cross 
section has the form [5l| : 



2v(V) 



dV 



cc 



dxdy 



TXS 



a 



2 sin 2 9 w M*y 



M 2 W 



! [xyFr(x,Q 2 ) 

,2 



+(i - y - x -^W\x, g 2 ) ± (y - y Y )xF.r(x, Q'% 



(28) 
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where Mz and M\y are the masses of the neutral- and charged- current interacting weak 
vector bosons, respectively, 9 W is the Weinberg angle, x = Q 2 /2p ■ q, y = p ■ q/p ■ k, and the 
square of the four momentum transfer for the reaction is Q 2 = —q 2 , k (p) is the momentum 
of initial state for neutrino or antineutrino (nucleon) and s = (k + p) 2 . Besides these, 



Z{W ± ) Vt 



x, Q 2 ) are the structure functions on proton (p), which only depend on x as Q 



oo, and are written in terms of the parton distributions as 

lim if p (x,Q 2 ) = l[((f-f + (fX) 2 )(u^x)+vF(x) + ( f(x)+^(x)) 

+ {Uvf + UlY) (d p (x) + 7{x) + s p (x) + s p (x) 
lim F 3 Zp (x,Q 2 ) = 2[f^(u p (x)-vF(x) + cf(x)-^(x)) 

Q — >oo 

+fvfl (d p {x) - d\x) + s p (x) - s p (x) 
F 2 Zp (x,Q 2 ) = 2xF Zp (x,Q 2 ). 

And the structure functions of charged- current have forms: 

lim F^ +p (x,Q 2 ) = d p (x) + u p (x) + s p (x) + c"(x), 



(29) 



lim FY p (x,Q' 2 ) = u p {x) + d P (x) + s p (x)+d'(x), 



lim if^(*,Q 5 



lim FY p (x,Q*) = u p {x)-d l '{x)-s p {x) + c p {x), 

2 Q 2 ^oo 



d p (x)-u p (x) + s p (x)-^(x), 



x,Q 2 ) = 2xFY p (x,Q 2 ). 



(30) 



One can obtain the structure functions for neutron n by replacing the superscripts p — m 
everywhere in Eqs. (J2HJ) and (j3~Uj) . In Eq. (J2HJ), fy, f%, fy and f\ are vector and axial- vector 
couplings: 

^ = I-l s in 2 ^, f\ = \ 



rd 

Jv 



- + - sin 2 6^, fi = — . 
2 3 2 



When the charge symmetry 



d n (x) 
u n {x) 
s n {x) 
c n {x) 



u p (x), 
d p (x), 

s p (x) = s(x), 
c p (x) = c(x), 



(31) 
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and c(x) = c(x) are valid, the structure functions for the isoscalar target with the s(x) ^ s(x) 
are given by 

lim F^ +No (x,Q 2 ) = hd p (x)+d P (x)+u p (x)+u p (x) + 2s(x) + 2c(x)}, 

Q 2 — >oo 2 

lim F^~ No (x,Q 2 ) = l[d p (x)+d P (x)+u p (x)+u p (x) + 2s(x) + 2c(x)}, 

lim F™ +N °(x, Q 2 ) = d p {x) + u p {x) - u p (x) - d\x) + 2s(x) - 2c(x), 

Q 2 ^oo 



lim F™ N °(x, Q 2 ) = d p (x) + u p (x) - u p {x) - d\x) - 2s{x) + 2c{x), 

Q 2 ~^oo 



°(x,Q 2 ) = 2xF™ JV °(x,Q 2 ), (32) 

where F™ +N °{x,Q 2 ) = \{F™ +P (x, Q 2 ) + F t w+n (x, Q 2 )). From above equations, we obtain 
the modified P-W relation with the asymmetry of six) and s(x) considered: 

a vN _ a VN c 

R N = „N _ „VN = R~ - $ R 7i (33) 

a cc a cc 

where 5R~ is the correction brought by the asymmetric distribution of strange-antistrange 
sea and R~ is the naive P-W relation with 

'K'V-WUchTW' (34) 

and we define Qy = J Q x[uv(x) + dy(x)]dx and S~ = x(s(x) — s(x))dx. From above 
equations, we can find that the asymmetry of s(x) and s(x) distributions should bring 
correction to sin 2 6 W . The results are given in table ITT1 and IIHI with the values of 5R~ ranging 
from 0.00297 (0.00312) to 0.00826 (0.00867) corresponding to A K = A n = 900 - 1500 MeV 
for different inputting parameters. However, the SUf(3) symmetry breaking requires smaller 
(Pk) and (P v ) j^, so that we should adopt a smaller value for such as from 900 MeV 
to 1100 MeV. Thus, the correction brought by the asymmetry of s(x) and s(x) distributions 
should remove the NuTeV anomaly by about 60-100% within the effective chiral quark 
model. Although the magnitude of S~ with CQ model as input is almost twice larger than 
that with CTEQ parametrization as input, as pointed out in Sec. m the calculated results 
of 5Rj given in Table HP and IIHI are similar and not sensitive to the different inputs at fixed 
A#. The reason may be that the uncertainties in the numerator and denominator of the 
Eq. (JHll) brought by the model can cancel out each other. 
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TABLE II: The results for CQ model input 





Z 


Qv 


s- 


5R- 


900 


0.74888 


0.86376 


0.00558 


0.00297 


1000 


0.73996 


0.85484 


0.007183 


0.00384 


1100 


0.73063 


0.84551 


0.00879 


0.00473 


1200 


0.72107 


0.83595 


0.01040 


0.00562 


1300 


0.71143 


0.82632 


0.01198 


0.00651 


1400 


0.70181 


0.81669 


0.01353 


0.00740 


1500 


0.69227 


0.80715 


0.01503 


0.00826 



TABLE III: The results for CTEQ parametrization input 





Z 


Qv 


s- 


SR 7 


900 


0.74888 


0.37089 


0.00252 


0.00312 


1000 


0.73996 


0.36686 


0.00322 


0.00402 


1100 


0.73063 


0.36247 


0.00398 


0.00498 


1200 


0.72107 


0.35831 


0.00468 


0.00590 


1300 


0.71143 


0.35395 


0.00539 


0.00683 


1400 


0.70181 


0.34960 


0.00612 


0.00780 


1500 


0.69227 


0.34528 


0.00675 


0.00867 



IV. DISCUSSION ON ADDITIONAL SYMMETRIC SEA CONTRIBUTIONS 

The asymmetries of light-flavor antiquark and strange-antistrange distributions originated 
from the nonperturbative sea have been discussed within the effective chiral quark model. 
Although the asymmetries of light-flavor antiquark and strange sea content in the nucleon 
are mainly coming from the intrinsic sea, the extrinsic sea should be also considered from a 
strict sense when we investigate the distributions of d(x) — u(x), d(x)/u(x) and s(x)/~s(x). 
From the calculations, we find that the effective chiral quark model predictions of d(x) —u(x) 
distributions keep consistence with the experimental data, but the predictions of d(x) /u(x) 
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distributions are not good in matching the experimental data at small and large x by us- 
ing the effective chiral quark model before taking into account additional contribution from 
the symmetric sea (e.g. 5d(x) = 5u(x)). Because the contribution of the symmetric sea 
effect (5d(x) = 5u(x)) cancels out each other in the distributions of d(x) —u(x), the distri- 
butions of d(x) —u(x) match the experimental data well. For the distributions of d(x)/u(x), 
when the symmetric sea content is considered, the distribution of d(x) /u(x) becomes: 

d(x) d(x) + Sd(x) — u(x) — 5u(x) 
u(x) u(x) + Su(x) 

= ^ ~w \ + !" ( 35 ) 

From the Eq. ([35)1 . we can know that the ratio of d(x) /u(x) will decrease if the symmetric sea 
is not zero. Here, we estimate the symmetric sea contribution 5u(x) + 5d(x) (5u(x) = 5d(x)) 
by the difference between the parametrization j3| of u(x) + d(x) and the results predicted 
by effective chiral quark model. Thus, the distributions of d(x)/u(x) without and with 
symmetric sea contribution are obtained and shown in Fig. |H1 From the figure, we can find 
that a large suppression occurs by including the symmetric sea contribution, which implies 
that the effective chiral quark model only provides a fraction of the total light-flavor in the 
nucleon, and that a significant fraction of symmetric sea content is required. This also can be 



found in Fig. [3 which disp 
from the parametrization 



ays the comparison of the distributions for u(x) and d(x) obtained 
and the effective chiral quark model calculations. The role of 
the symmetric sea is similar to that by including some isoscalar (pseudoscalar) meson a (rj) 
contribution in the meson cloud model, as pointed in Ref. 5^|. However, the effective way 
of taking into account the symmetric sea contribution is not good in reproducing the shape 
of d(x)/u(x), and more work both theoretically and experimentally should be done along 
this direction. 

Within the effective chiral quark model, the strange-antistrange asymmetry mainly comes 
from the K meson clouds which accompany the constituent quarks. From the calculations 
we know that the corrections coming from the strange asymmetry with different inputs of 
constituent quark distributions are large enough to explain the NuTeV anomaly. It is also 
necessary to mention that the ratio of s(x)/~s(x) in the model does not necessarily conflict 
with the experimental data in fact, because the experimental data include the effect of the 
intrinsic and extrinsic sea. In this paper, we calculate the ratio of s(x)/s(x) without and 
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FIG. 6: The distributions of d(x)/u(x) for A n = 1500 MeV. The thick solid and dotted curves are 
the chiral quark model results with CTEQ6 parametrization and the CQ model as inputs. The thin 
curves are the corresponding results with both chiral quark model results and the symmetric sea 
contribution Su(x) = 5d(x) estimated by difference between the NuTeV data parametrization and 
the chiral quark model results. And the experimental data come from E866/NuSea (Q 2 = 54 GeV 2 " 1 
experiments 



UK i 

3- 



with symmetric strange sea contribution, where the symmetric sea contribution is estimated 
effectively by the difference between s(x) +s(x) from the NuTeV data parametrization and 
that in the effective chiral quark model calculation, in similar to what we estimated above 
for the symmetric light-flavor sea quarks. The so called "NuTeV data parametrization" here 
is obtained by the NuTeV collaboration from the results of the leading order fits to the cross- 
section extracted from the NuTeV experimental data . The strange sea parameterizations 
which break s(x)dx = s(x)dx as pointed out as in Q, 1^ are defined conventionally 
by 

u{x) + d{x) 
s(x)=k (1-aO , 

. . u(x) + d(x) ,„ .- . . 

s(x) = k K ' 2 1 \ l-x) a . (36) 

where the values of k, k, a, a are taken from the NuTeV set of the Table III in Ref. [5^. And 
we give the comparison of s(x)/s(x) from the model prediction and that from the dimuon 
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FIG. 7: The distributions of d(x) and u(x) for A,f = 1500 MeV. The solid and dotted curves are 
the effective chiral quark model results with CTEQ6 parametrization and CQ model inputs, and 
the dashed curve is for the NuTeV data parametrization. 



measurements JscJ in Fig. |HJ From the left one of the figure without the symmetric sea 
content, we can see that the ratios are out of the area for errors. However, the distributions 
of s(x)/~s(x) including the symmetric sea contribution are almost within the range of the 
experimental errors in the right one. This case is similar to the d(x) /u(x) distribution and 
indicates that the effective chiral quark model only provides a fraction of the strange sea 
content, and that a significant fraction of symmetric strange sea content is still needed. 
Fortunately, the strange asymmetry of the effective chiral quark model can cause a large 
reduction of the NuTeV anomaly, although the asymmetric strange sea is only a small part 
of the total strange sea content of the nucleon. 
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FIG. 8: Distributions of s(x)/s(x), where the shadowing area is the error range of NuTeV. The 
thick and thin solid curves are the effective chiral quark model results for s(x)/s(x) when Ak = 900 
and 1100 MeV with CTEQ6 parametrization as input. The thick and thin dotted curves are the 
effective chiral quark model results for s(x)/s(x) when Ak = 900 and 1100 MeV with CQ model 
as input. The left side is the prediction by the effective chiral quark model only and the right side 
is the result by including both the prediction of the effective chiral quark model and the symmetric 
sea contribution estimated by the difference between the NuTeV data parametrization and the 
effective chiral quark model result. 

V. SUMMARY 

In this work, we presented the calculations of the asymmetries of light-flavor sea quark 
distributions d(x) — u(x) and strange-antistrange s(x) — s(x) within the effective chiral quark 
model with more details along with our previous work |22j by using two different sets of 
parameterizations as inputs and find that the results for d(x) —u(x) match the experimental 
measurements well. The distributions of d(x)/u(x) do not match with the experimental 
data, but there is a large suppression of the ratio when the symmetric sea contribution 
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is considered. The contributions of symmetric sea effects are estimated effectively by the 
difference between the effective chiral quark model results and the data parametrization 
results. Fig. 0] and Fig. [7| indicate that the predictions of the effective chiral quark model 
provide only a fraction of the total light-flavor sea quarks and strange sea content of the 
nucleon. Also, we point out that the calculated results for s(x)/s(x) in the effective chiral 
quark model are consistent with the available experiments with an additional symmetric 
sea contribution being included effectively. More noticeably, the asymmetry of strange- 
antistrange distributions can bring a significant contribution to the NuTeV defect by at 
least 60% with reasonable parameters within this model, and the results are not sensitive 
to the different inputs of constituent quark distributions, although the intrinsic sea is not 
the dominant part in the nucleon sea. The effective chiral quark model is thus successful in 
explaining the nucleon sea anomalies not only about the GSR violation and the proton spin 
problem, but also about the NuTeV defect. Also, this may imply that the NuTeV anomaly 
can be considered as a phenomenological support to the strange-antistrange asymmetry 
of the nucleon sea. So it is more important that more precision experiments should be 
carried out to enable direct and more accurate determination of strange-antistrange sea 
quark distributions in the future. 
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